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Abstract. We report resonant features of novel controllable reflectarray which 
consists of meander-like graphene strips placed on a metal-backed dielectric substrate. 
The structure manifests two kinds of resonances appeared as sharp deeps of reflectivity. 
The first one exists because the strips of periodic cells of the structure have resonant 
sizes for induced surface plasmon-polaritons. The second kind of resonances is defined 
by excitation of TM eigenwaves of the whole structure as a plane photonic crystal. The 
latter resonances do not depend on whether the strips of the unit cells have resonant 
sizes or not. 
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1. Introduction 

Planar arrays of conducting elements find 
many applications, especially as frequency se¬ 
lective surfaces (FSS) and slow-wave struc¬ 
tures. In recent years, periodical structures 
with different types of dielectric and conduct¬ 
ing elements became very popular in the design 
of the so-called electromagnetic crystals also 
known as photonic band gap structures p. 

Two-periodic planar strip arrays are 
promising for applications because they can 
possess resonant properties in frequency region 
of a single-wave regime due to complex shape 
of their elements. Besides, they have very 
small thickness. Transmission and reflection 
properties of the complex-shaped metal strip 
arrays were studied recently with reference to 
potential applications of them as hlters, thin- 
hlm sensors and slow-light devices (see, for 
example, a review paper 0)- In many of these 
works, it was supposed that the strip particles 
placed in different cells of a periodic structure 
have no electric connection. 

The purpose of this work is to study 
resonant reflection properties of 2D (two- 
dimensional) periodic array consisting of long 
continues graphene wavy-shaped strips [3] 
(“hsh-scale” HI), placed on a metal-backed 
dielectric substrate. A unit cell of the array is 
shown in Fig. [1^. We report resonant features 
of this array as a planar photonic structure 
which were not investigated before. 

The hsh-scale morphology, also called 
meander-line structure possesses resonant 
properties. Such metallic gratings are com¬ 
monly used in optical magnetic mirrors and 
polarization transformers m i- They can 
serve also as hlters [7]. There are ideas to ex¬ 
ploit electromagnetically induced transparency 
in such structures [8]. 

Graphene is a 2D material with only 


one atom thickness which manifests unique 
electronic and electromagnetic properties |9l 
[in]. Electrically tunable carrier concentration 
of graphene is used in hybrid devices, where 
metal plasmonic structures and graphene 
are combined to broadening and tuning 
resonances m- Besides, plasmonic resonances 
directly in the graphene elements provide 
a great potential for designing devices with 
a remarkably high absorption [la US] and 
unprecedented levels of held conhnement |14j . 
Several pure graphene planar arrays with 
diherent functionality have been suggested 
recently [I511I6]. 

When the graphene structure is placed 
on a metallic mirror coating with a thin 
dielectric layer as a substrate it becomes a good 
broadband rehector everywhere apart from 
few frequency regions where the rehectivity 
is small. At these frequencies the rehected 
wave shows no or a small phase change 
with respect to the incident wave, thus 
resembling a rehection from a hypothetical 
magnetic wall. We also observed that the 
structure acts as a local held concentrator 
and a resonant ampliher of losses in graphene 
strips and underlying dielectric layer. Control 
bias potential in the hsh-scale graphene 
structure discussed below can be applied 
via uninterrupted strips. The bias voltage 
changes chemical potential of graphene and 
consequently, its conductivity. 

We will show that there are two kinds of 
resonances in hsh-scale graphene rehectarray 
which diher by physical nature. The hrst 
one exists due to strips of the periodic 
cell having resonant dimensions for excited 
surface plasmon polaritons (SPPs). In the 
case of metal-made hsh-scale array, these 
kind of current resonances is known and well 
investigated |1]. 

The second kind of resonances of this 
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array results in excitation of eigenmodes of the 
whole structure as a planar photonic crystal. 
These resonances almost do not depend on 
whether the graphene strips of the unit cells 
have resonant sizes or not. To the best of our 
knowledge, the resonances of this kind in hsh- 
scale array were not investigated before. 

2. Description of fish-scale refiectarray 

Fish-scale graphene array is patterned on the 
interface of dielectric substrate which has 
the thickness h and non-dispersive relative 
permittivity £ = 3.5. A unit cell of the 
array is shown in Fig. [1^. The graphene strip 
has the width w, its thickness is inhnitely 
small in comparison with all structure sizes 
and wavelength. For the sake of simplicity, 
a square shape of unite cell is chosen. The 
pitch of the array is I = lOO/im. A perfect 
electric conductor (PEC) screen is placed on 
the opposite side of the substrate. 
and k are electric held, magnetic held and wave 
vector of the incident plane wave, respectively. 

3. SPPs supported by graphene placed 
on PEC-backed dielectric substrate 

To study and explain the physical nature of 
resonant features of graphene rehecting hsh- 
scale array, we start with a discussion of basic 
properties of plasmon-polaritons propagating 
in investigated layered structure without 
patterning of graphene sheet. 

In problems of macroscopic electromag¬ 
netics in low THz frequency region, the elec¬ 
tric properties of graphene can be described 
by surface conductivity m 

_ e^ksT r He 
TTh?‘{u: — i2T) ksT 
where /Xc is the chemical potential, T is 
the temperature (the value of 300 K was 
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Figure 1. (a) Unit cell of fish-scale graphene array 
placed on PEC-backed dielectric substrate, coordinate 
frame, and main sizes: I = 100 /xm, w = 30 ym, h = 
25 ^m, yc = 0.5 eV, resonant frequency is 0.45 THz. 
(b) Surface current distribution in case of excitation 
of resonance of dimensional kind by normally incident 
y-polarized wave. 
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Figure 2 . Surface conductivity of graphene Recr 
(red solid line) and Imtj (green dashed line) versus 
frequency, Uc = 0.5 eV. 

chosen in the simulations), F = ( 2 r)“^ is the 
scattering rate, r is the relaxation time (set 
as 1 ps). Frequency dependence of graphene 
conductivity is shown in Fig. |2l 

Let us consider TM polarized SPP 
supported by graphene. A basic structure is 
illustrated in the insert of Fig. HJ According to 
the classical waveguide theory, magnetic held 
of the SPP propagated along a:-axis in different 
regions of the layered structure can be written 
as follows: 

Hi = z>0 

H2 ^ + C'e-*^ 2 ( 2 +fe)]^ -h < z < 0 , 

where 71 = \/k‘^ — and 72 = \/k‘^e — 
while Im 7 i < 0 and Im 72 < 0 with fi being 
propagation constant of the SPP, k is the wave 
number in free space. By using the boundary 
conditions, that is, the transverse to 2 ;-axis 
component of electric held is zero in the plane 
z = —h and continuous in the plane z = 0 
but the transverse to z-axis magnetic held has 
a discontinuity in this plane which is equal 


to the value of surface current density on the 
graphene sheet: 

Hl-Hl = -aE,, z = 0, (3) 

the dispersion equation of SPP can be derived: 
7i£(l+e-“’“'‘)+72(l+<T.)|-)(l-e-“^"'‘) = 0, (4) 

where rj is the impedance of free space. 

In case of the space h between the 
graphene sheet and PEC is inhnitely large, the 
dispersion equation for the SPP propagating 
along graphene placed on a dielectric semi¬ 
space follows from (II]) and can be written as 

7i£ + 72(1-h cr?7^) = 0 . ( 5 ) 

The well known expression for the propagation 
constant for the SPP supported by a free 
standing graphene layer [HI follows from 
equation ([5]): 

z? = (6) 

where we put e = 1. The frequency depen¬ 
dences of the corresponding SPP propagation 
constants are shown in Fig. [3] and Fig. HJ 

Lastly, in case of the conductivity a equals 
to zero, i. e. the interface = 0 is free from 
graphene, the dispersion equation for the TM 
eigenwaves follows from ( 01 ) and can be written 
as 

\/k‘^e — /3^tan {h^Jk?e — f3‘^) = e \/— k^.ifi) 

PEC-backed dielectric layer is non dissi¬ 
pative, at least one real solution of equation 
d?]) exists corresponding the fundamental low¬ 
est frequency TM mode which has no cut-off 
frequency. The frequency dependence of the 
mode propagation constant is shown in Fig. [Sj 

4. Analysis of resonant properties of 
graphene refiectarray 

In this Section, we consider a case of array 
excitation by normally incident ^/-polarized 
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Figure 3. Frequency dependences of real parts of 
propagation constants of SPPs in cases of free standing 
graphene sheet (dotted line), graphene placed on 
interface of dielectric semi-space which permittivity is 
e = 3.5 (red dashed line) and on surface of PEC-backed 
dielectric substrate for two different thickness 25 ym 
(blue solid line) and 10 ym (green dot-dashed line), 
yc = 0.5 eV. 
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Figure 4. Frequency dependences of imaginary parts 
of propagation constants of SPPs in same cases as in 

Fig.H 



Figure 5. Propagation constant of fundamental TM 
mode of PEC-backed dielectric layer (e = 3.5) of 25 ym 
thickness versus frequency. 

plane wave (see Fig. [1^). In the conclnsions 
we shall give short comments related to x- 
polarization of the wave. From the dispersion 
dependences presented in Fignres[3l [Hone can 
see that retardation of the eigenwave of non- 
patterned graphene (on-dielectric PEC-backed 
strnctnre) is larger for the dielectric layer with 
smaller thickness. Bnt redncing the strnctnre 
thickness results in increasing of the wave 
damping (see Fig. 0]). Thus we have chosen 
the dielectric substrate thickness h = 25 pm 
which is typical for THz devices. 

4-1. Features resulted from resonant 
dimensions of graphene strips 

The reflectivity of structure manifests in a 
resonant deep resulted from excitation of SPPs 
with wavelength which is approximately equal 
to the length of graphene strip inside a unit 
cell. As a consequence of the structure 
symmetry, the incident y-polarized wave does 
not excite currents at the edge points of the 
strip where the electric field is orthogonal to 
the strip. To estimate the value of resonant 
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frequency, we compare the length of the line 
drowned between points A and B (they are 
placed on the boundaries of unit cell, see 
Fig. [1^) with zero values of surface current 
density and the wavelength of SPP by data of 
solution of dispersion equation dH) (see Fig. [3]). 
Approximate resonant condition is reduced to 
equality of the distance between points A and 
S to a half of the SPP wavelength. This 
estimation leads to the resonant frequency 
0.563 THz. 

Frequency dependence of reflection coeffi¬ 
cient obtained by full waves numerical simula¬ 
tions is presented in Fig. [6l If one takes into 
account a complexity of the shape of the ar¬ 
ray elements and the essential decay factor of 
SPPs, one can see that our analytical estima¬ 
tion of the resonant frequency and the related 
numerical data are in a good agreement. Thus, 
this lowest frequency resonance of the reflec- 
tarray depends on the strip dimensions and it 
is similar to resonances of classical microstrip 
resonators and antennas in microwave region 
[18] and also resembles resonances of longitudi¬ 
nal currents in the hsh-scale array with narrow 
PEC strips [5]. In our case, the resonant sur¬ 
face current distribution is shown in Fig. [TJj. 

The propagation constant of SPPs de¬ 
pends on graphene chemical potential, there¬ 
fore the resonant frequency of the structure can 
be changed by an external bias voltage. In the 
insert of Fig. [6] we present the dependence of 
the resonant frequency on the chemical poten¬ 
tial Pc- 

4-2. Resonances resulted from excitation of 
eigenmodes 

In periodic structures, the scattered held of 
normally incident plane wave has discrete 
space spectrum and can be presented as a su¬ 
perposition of the principal waves propagating 



Figure 6. Frequency dependence of reflection 
coefficient for h = 25 gm, ru = 30 ^m, 

Uc = 0.5 eV, normally incident j/-polarised wave. 
Insert: dependence of resonant frequency on chemical 
potential gc- 

along ; 2 -axis and partial waves (propagating 
or evanescent) both inside the substrate and 
above the graphene in the near free space. The 
lowest partial waves propagate in plane xOz 
(one pare of waves) and yOz (another pare). 
They have identical absolute value of x- and 
y-components of the wave vectors 271/I. 

To observe resonances of rehection based 
on energy store in electromagnetic held 
oscillation resulted from excitation of opposite 
propagating eigenwaves, it is necessary to fulhll 
the following requirements. In the plane xOy 
the propagation constants of both eigenwave 
of the structure (3 and the partial wave of the 
scattering held must be identical: 

fiif) = 277/1. (8) 

It is a synchronism condition. 

Next, to have a high quality factor of 
resonance, it is necessary to exclude energy 
loss via in plane propagating eigenwaves. Thus 
the amplitudes of excited eigenwaves along 
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some direction and in opposite one must be 
identical. Let us consider firstly the excitation 
of eigenmodes propagating in a:-direction. The 
amplitudes of excited eigenwaves along Ox 
axis and in the opposite direction must be 
equal. This requirement is fulhlled in case of 
a hsh-scale structure because the plane yOz 
is its symmetry plane. In case of eigenwaves 
propagated along ^/-direction, the condition is 
more complex because the plane xOz is a glide 
symmetry plane of the array. However, it 
is sufficient condition to exclude a direct or 
reversal power flow. 

Thirdly, the phases of the opposite 
propagating eigenwaves must correspond to 
constructive interference of their helds inside 
the substrate with the scattered held of the 
incident wave. This condition is satished also 
due to the shape of graphene strip in the unit 
cell, as it can be seen easily. For instance, let 
us consider the case of eigenwaves propagated 
along a:-axis. The longitudinal current has zero 
values in the points dividing strip inside a unit 
cell in two parts (they are the points A and 
B shown in the Fig. [1^). These two parts 
are F-shaped planar chiral particle and mirror 
to S'-shaped one with opposite handedness. 
Thus, on the surface of these particles the 
longitudinal currents induced by incident y- 
polarized held excite magnetic helds having 
opposite directed Hy. These helds are shifted 
in //2 along x-axis that is the same as right- 
and left-handed S'-shaped particles inside unit 
cell (see Fig. [7|i). 

It is also worth to mention that the 
main held components of the excited resonant 
oscillation which are Hy^ and E^, are absent 
in the incident plane wave. As a consequence, 
coupling of the oscillation and the incident 
wave is small and the radiation loss is also 
small. Thus we deal with the trapped mode 
regime of the open planar periodic structure. 


E inc 



(b) 


Figure 7. a) i/y-component of resonant oscillation 
(/=1.887 THz), eigenwaves propagate in x-direction, 
excitation by y-polarized normally incident wave, h = 
25 /rm, w = 15 ym, /ic = 0.5 eV. b) Surface current 
distribuition in this case 

For the sake of brevity we will not now consider 
in detail the case of resonant excitation of eigen 
waves propagated in ^/-direction. Thus the hsh- 
scale patterned structure is a perfect object 
to observe resonant features resulted from 
excitation of eigenmodes of planar photonic 
crystal. 

The spectrum of rehection has been 
obtained numerically and it is shown in 
Figures [HI [H As one can see, the resonances 
of these kind are excited by ?/-polarized wave 
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Figure 8. Frequency dependence of reflection 
coefficient: h = 25 /rm, fXc = 0.5 eV, normally incident 
2 /-polarized wave. 



(a) 




Figure 9. Frequency dependence of phase of reflection 
coefficient: h = 25 ym, = 0.5 eV, normally incident 
y-polarized wave. 


Figure 10. i/j^-component of resonant oscillation, 
/=2.031 THz, eigenwaves propagate in y-direction, 
excitation by y-polarized normally incident wave, h = 
25 /im, re = 15 ym, yc = 0.5 eV. b) Surface current 
distribuition in this case 

as a doublet for every of the chosen strip width. 
The resonant frequencies and the difference 
between them are collected in Table 1. The 
shown data reveal that the narrow strips have 
close resonant frequencies in doublet because 
the features of the array with narrow strips 
are closer to the properties of PEC-backed 
dielectric as a waveguide. 

The lower resonant frequency of each 
doublet corresponds to the excitated eigen¬ 
waves propagating along x-direction and the 
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Table 1. Resonant frequencies of doublet for different Table 2. Resonant frequencies and absorption of 
w, e = 3.5, h=25 ym, yc = 0.5 eV doublet for different yc, e = 3.5, h = 25 yui, w = 

15 ym 


w 

/i |THz] 

/2 [THz] 

A/ = /2 - /i [THz^ 







5 

2.049 

2.076 

0.0274 

/ic [eV] 

/i |THz| 

/2 [THz] 

Ai 

[%] 

A 2 

[%] 

10 

1.986 

2.070 

0.084 

0.3 

1.722 

2.142 

99.6 

70.8 

15 

1.887 

2.031 

0.144 

0.5 

1.887 

2.031 

98.2 

99.3 





0.8 

1.923 

2.118 

98.3 

99.3 





1.0 

1.980 

2.166 

65.6 

84.9 


higher frequency corresponds to propagating 
ones along t/-axis. As the main component of 
plasmon-polariton magnetic field at lower reso¬ 
nant frequency is orthogonal to the same com¬ 
ponent of the incident wave, the quality factor 
in this case is higher than the one related to 
the second resonance where these components 
are parallel resulting in a bigger radiation loss. 

The transverse to strip surface current 
components prevail at the excitation of this 
kind of resonances (see Figures [7 )d andfTOb ) in 
contrast to the longitudinal currents at the res¬ 
onances considered above in section 4.A. Be¬ 
sides, in the second resonance of doublet, be¬ 
cause the structure has glide symmetry plane 
being orthogonal to propagation direction of 
the eigenwaves, the related standing wave 
manifests a complex held distribution inside 
the unite cell. It is resulted in appearance of 
regions on graphene strips with specihc rotat¬ 
ing current distributions which are shown in 
Fig. [TUb . 

The quality factor of resonance depends 
on dissipative properties of the materials of the 
structure and on diffraction loss. For the sake 
of simplicity, we assumed that the dielectric 
substrate is nondissipative. 

Features of eigenwave of the considered 
structure differs in a small way from those 
of eigenwave of PEC-backed substrate in the 
case of narrow graphene strips. The resonant 
frequencies are far from those of the above 


studied dimensional resonances of the strips. 
Thus, the frequency of this photonic crystal 
kind of resonances can be estimated by using 
in the equation (|8]) instead of rigorous value 
/5, an approximated one related to eigenwave 
of backgrounded substrate and presented in 
Fig. |5l The solution of approximate 
equation ([8]) provides estimation of resonant 
frequency as / = 2.057 THz which is in 
good correspondence with resonant frequencies 
following from numerical simulations (see 
Table 1). The frequency discrepancy is smaller 
when the strip width is narrower because 
perturbation of the eigenwaves is smaller. 

In Table 2 we present data characterizing 
dependence of the resonant frequencies of the 
doublet and absorption (A = 1 — |5'iip) 
on the chemical potential fic- One can 
see a smaller variation of resonant frequency 
in comparison with the discussed above 
dimensional resonance tuning. It is evident 
consequence of quite different physical nature 
of these kinds of resonances. However, by 
changing chemical potential due to variation 
of bias voltage, one can modify the level of 
absorption of the structure. It is important for 
application in devices for modulaton of signals 

PI- 
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5. Conclusions 

We reported features of novel tunable graphene 
fish-scale reflectarray. Uninterapted strip 

elements of array are convenient to apply 
control bias voltage. The structure manifests 
two kinds of resonances which differ by 
physical nature. Analysis reveals that they 
are excited by incident wave in cases of 
polarizations both along the strips of array 
and transverse to them. The first one 
exists because inside of periodic cells the 
strips of structure have resonant lengths 
for induced currents. The reflectivity of 
structure manifests in a resonant deep (peak 
of absorption) resulted from excitation of SPPs 
with wavelength which is approximately equal 
to the length of graphene strip inside a unit 
cell in the case of excitation by incident wave 
polarized transversally to array strips. 

The second kind of resonances of fish-scale 
array placed on the metal-backed substrate is 
resulted from excitation of TM eigenwaves of 
the structure which have identical amplitudes 
but propagate in opposite directions. An 
energy is stored in the standing wave along 
the whole structure being a plane photonic 
crystal. These resonances are not dependent 
on whether the strips have resonant sizes 
within the unit cell or not. 
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